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SUMMARY 


The  work  summarized  in  this  report  covers  the  first  quarter  of  a  program  with  a  goal  that 
is  twofold:  first,  to  study  the  properties  of  native  point  and  line  defects  in  IRFPA  active  and 
substrate  materials,  and  second,  to  study  the  properties  of  native  point  defects  in  two  classes  of 
photonic  materials,  the  wide-gap  II-VI  compound  (ZnSe  as  the  prototype  for  which  impurity 
properties  will  also  be  calculated)  and  the  nonlinear  optical  materials  (LiNbOs  as  the  prototype). 
Our  accomplishments  in  the  first  quarter  include: 

•  Preliminary  determination  of  the  concentration  of  native  point  defects  in 
Hgo.8Cdo.2Te  as  a  function  of  anneal  temperature  and  mercury  pressure. 

•  Identification  of  mechanisms  likely  to  account  for  effects  of  dislocations  on  IRFPA 
performance,  and  initial  quantitative  analysis. 

•  Identification  of  features  of  the  band  structure  important  to  native  point  defect 
determination  and  also  important  to  transport  and  lifetime  analyses. 

•  Initial  calculation  of  the  properties  of  ZnSe,  including  new  gradient  correction  terms. 
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1. 


INTRODUCTION 


We  have  added  a  number  of  improvements  to  our  calculation  of  the  native  defect  forma¬ 
tion  energies;  these  improvements  will  apply  to  HgCdTe,  ZnSe,  and  LiNbOs.  First,  we  have 
included  a  gradient  correction  to  the  local  density  approximation  that  impacts  the  manner  in 
which  the  exchange-correlation  energy  is  calculated.  This  represents  a  signiHcant  correction  to 
energies  in  which  atoms  are  exchanged  between  the  solid  and  free-atoms  state;  for  example,  with 
this  correction  the  cohesive  energies  of  the  II-VI  compounds  are  predicted  in  excellent  agreement 
with  experiment  (<  5%). 

We  have  also  begun  to  calculate  the  full  nonradial  relaxation  of  the  defects;  because 
relaxations  may  be  large,  some  interatomic  distances  increase,  while  others  decrease.  This  has 
necessitated  careful  examination  of  parameters  used  in  the  codes,  such  as  the  size  of  the  basis  set. 
We  are  completing  this  analysis  and  have  begun  to  calculate  relaxation  energies. 

Finally,  we  have  developed  an  approach  that  has  permitted  us  to  identify  the  primary  ion¬ 
ization  state(s)  of  defects.  Although  the  local  density  approximation,  upon  which  our  methods 
are  based,  predicts  band  gaps  that  are  too  small  in  the  semiconductors,  this  approach  appeals  to 
yield  reasonable  agreement  with  the  available  experiments  for  the  prototypical  defect  we  have 
studied.  Later  in  the  program,  we  will  be  making  improvements  to  this  portion  of  our 
calculations. 


2.  NATIVE  POINT  DEFECTS  IN  IRFPA  ACTIVE  AND 
SUBSTRATE  MATERIALS 


2.1  NATIVE  POINT  DEFECTS  IN  Hgo^CdO^Te 

We  have  included  the  gradient  correction  in  the  native  defect  formation  energies  for 
HgTe,  and  calculated  the  resulting  defect  densities  in  Hgo.8Cdo.2Te  as  a  function  of  temperature 
and  mercury  pressure.  The  ionization  states  of  the  defects  have  also  been  included  in  our 
analysis.  We  find  good  agreement  with  the  available  experimental  data. 

We  confirm  the  mercury  vacancy  as  the  dominant  defect  in  narrow-gap  HgCdTe.  We 
predict  that  the  tellurium  antisite  will  be  an  important  secondary  defect,  and  are  currently  inves¬ 
tigating  if  it  is  the  “residual  donor”  in  annealed  HgCdTe.  The  results  of  our  work  were  presented 
at  the  1992  US.  Workshop  on  the  Physics  and  Chemistry  of  Mercury  Cadmium  Telluride  and 
Other  IR  Materials,  held  in  Boston  on  13  through  15  October.  A  copy  of  the  paper  presented  at 
that  meeting  is  presented  in  Appendix  A. 
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Because  HgCdTe  is  grown  and  processed  at  temperatures  at  which  the  material  is  intrin¬ 
sic  and  degenerate,  and  because  the  defects  are  ionized,  too,  and  interact  with  the  electron-hole 
gas,  we  must  know  the  intrinsic  reaction  constant  in  order  to  calculate  the  concentration  of  ion¬ 
ized  defects.  While  there  has  been  much  study  of  the  band  edge  properties  of  HgCdTe  at  low 
temperatures,  the  relevant  properties,  (band-gap  Eg,  and  electron  and  hole  effective  masses)  at 
the  elevated  temperatures  have  not  been  studied.  In  the  paper  in  Appendix  A,  we  used  expres¬ 
sions  for  Eg  fix>m  the  literature,  and  extrapolated  them  to  the  high  temperatures;  temperature 
independence  of  the  shape  of  the  conduction  and  valence  bands  was  assumed.  Such  an  extrapo¬ 
lation  must  be  replaced  by  better  experimental  results  or  further  theory.  In  addition,  we  have 
found  that  the  Moss-Bomstein  shift  is  large  at  the  higher  temperatures  (of  the  order  of  several 
tenths  of  an  eV),  and  needs  to  be  corrected  for  in  the  absorption  results  in  the  literature. 

Furthermore,  we  have  found  from  our  theoretical  calculations  that  the  electron  band 
shape  is  highly  nonparabolic,  except  very  near  to  the  band  edges.  The  shape  of  the  bands,  and 
their  temperature  dependence,  must  be  known  in  order  to  calculate  the  intrinsic  reaction  constant 
at  high  temperatures.  Even  low-temperature  properties  for  n-type  material  with  n  >  10'^  will 
have  interpretations  much  different  from  those  deduced  from  parabolic  band  structures.  These 
modified  properties  include  the  electron  mobilities,  lifetimes,  absorption  coefficients,  responsivi- 
ties,  and  noise  currents. 

2.2  DISLOCATIONS  IN  HgCdTe 

A  large  body  of  experimental  evidence  has  established  that  dislocations  adversely  affect 
the  device  performance  of  photovoltaic  detectors.  It  is  generally  agreed  that  the  presence  of  dis¬ 
locations  r^uces  the  minority  carrier  lifetime ^*2  and  diffusion  length^*^  of  the  semiconductor, 
increases  the  leakage  current^,  and  reduces  the  resistance^  and  the  open  circuit  voltage^  of  the 
detector  diodes.  Johnson  et  al.^  have  very  recently  summarized  the  effect  of  dislocations  on  the 
electrical  and  optical  properties  of  HgCdTe  (MCT),  and  Shin  et  al.^  have  explored  techniques 
for  reducing  dislocation  densities  in  MCT  grown  heteroepitaxially  on  GaAs  and  Si  using  minor¬ 
ity  carrier  density  as  a  diagnostic  tool. 

Dislocations  are  thought  to  affect  the  properties  of  MCT  by  acting  as  scattering  and 
recombination  centers.  To  account  for  this  activity,  it  is  necessary  to  explain  how  dislocations 
can  give  rise  to  the  electric  fields  that  affect  carrier  dynamics.  Such  electric  fields  can  be  asso¬ 
ciated  with  the  band  bending  arising  from  gap  states  in  dislocation  cores,  and  from  the  strain 
frelds  in  the  dislocation’s  vicinity.  We  are  considering  both  these  effects  for  the  present  study, 
because  it  appears  that  both  may  be  needed  to  account  for  the  observations,  and  for  the  develop¬ 
ment  of  effective  design  rules  for  MCT  photovoltaic  detector  arrays. 

Dislocation  cores  must  certainly  give  rise  to  bound  states  in  the  semiconductor  band  gap, 
and  it  is  believed  that  those  dislocation  cores  terminating  on  a  line  of  Te  atoms  in  n-type  MCT 
may  be  inverted,  and  surrounded  by  a  depleted  space-charge  region.  Analysis  of  the  fields  in  this 
region  using  the  depletion  approximation  predicts  depletion  region  radii  of  the  order  of,  or  even 
smaller  than,  the  Debye  length.  A  more  accurate  analysis  using  a  somewhat  less  crude  approxi¬ 
mation  (see  Appendix  B)  predicts  that  small  but  still  noticeable  fields  persist  to  larger  radii, 
resulting,  for  example,  in  a  radius  of  the  field  region  of  a  dislocation  in  n  =  10^^  cm-3  MCT  of 
the  order  of  5  ^im.  Minority  carriers  within  this  radius  will  have  their  lifetimes  shonened.  While 
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this  can  give  rise  to  significant  effects,  it  does  not  appear  adequate  to  account  for  the  observed 
phenomena. 

The  effect  of  strain  on  the  current- voltage  relation  of  InSb  and  MCT  has  been  investi- 
gated^.S  and  the  observations  have  been  interpreted  in  terms  of  the  piezoelectric  fields  associated 
with  the  strain.  Strain  scattering  by  point  defects  in  piezoelectric  semiconductors  has  been  ana¬ 
lyzed  by  Fedders.^  We  are  carrying  out  a  similar  analysis  (see  Appendix  C)  for  piezoelectric- 
induced  effects  caused  by  the  strain  field  of  dislocations.  Some  aspects  of  this  interaction  can  be 
understood  from  symmetry  considerations.  The  piezoelectric  tensor  of  the  sphalerite  structure 
couples  only  shear  strains  to  the  crystal’s  electric  polarization.  The  strain  field  of  a  long  straight 
dislocation  and  the  electric  polarization  that  it  induces  are  planar,  that  is  to  say,  independent  of 
the  coordinate  along  the  dislocation  direction.  The  piezoelectric  field  that  scatters  the  carriers  is 
generated  by  the  bound  charge  density  associated  with  the  divergence  of  the  polarization.  The 
polarization  has  nonzero  divergence  only  if  the  strain  displacement  (and  Burger’s  vector)  has  a 
component  along  the  dislocation  axis.  Therefore,  straight  edge  dislocations  produce  no  scatter¬ 
ing  field,  but  screw  and  mixed  dislocations  do.  Kinks,  jogs,  and  nonparallel  dislocations  are  also 
expected  to  generate  piezoelectric  fields  that  can  act  as  scattering  and  recombination  centers.  For 
a  quantitative  analysis  we  need  the  piezoelectric  coupling  constant'®  of  MCT  and  structural 
information  about  dislocation  types.  This  work  is  in  progress. 


3.  PHOTONIC  MATERIALS 


For  research  on  the  wide-gap  II- VI  compound,  ZnSe  was  chosen  as  the  prototype  and 
much  of  the  effort  has  been  in  establishing  the  groundwork  for  doing  the  relaxed  defect  calcula¬ 
tions  and  the  gradient  correction,  as  discussed  in  the  introduction.  The  cohesive  energy  and  bond 
length  of  ZnSe  have  been  calculated  and  are  found  in  good  agreement  with  experiment.  We  have 
begun  the  calculation  of  the  native  defect  formation  energies.  (Nation  and  anion  vacancies,  anti¬ 
sites,  and  two  classes  of  cation  and  anion  interstitials  are  being  included. 


4.  WORK  PLANNED 


In  the  next  quarter,  we  will  continue  work  to  improve  the  value  of  the  intrinsic  reaction 
constant  that  is  necessary  to  calculate  the  ionized  defect  concentration  in  HgCTdTe.  We  will  cal¬ 
culate  the  relaxed  neutral  native  defect  energies  in  HgCdTe  and  in  ZnSe.  We  will  examine  eight 
defects:  the  cation  and  anion  vacancies,  the  cation  and  anion  antisites,  and  the  cation  and  anion 
at  two  different  interstitial  sites.  Upon  completion  of  these  calculations  we  will  redo  the  ionized 
defect  density  calculations  in  HgCdTe,  and  make  our  first  predictions  of  neutral  native  defect 
densities  in  ZnSe. 
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Hgo.tCdo.iTe  Native  Defects:  Densities  and  Dopant  Properties 

M.A.  Berding,  M.  van  Schilfgaarde,  and  A.  Sher 
SRI  International,  Menlo  Park,  CA  94025 

We  examine  the  native  defect  equilibrium  in  HgCdTe,  including  cation  and  anion 
vacancies,  interstitials,  and  antisites  in  the  analysis.  A  gradient  correction  to  the  local 
density  functional  has  been  added  to  the  defect  formation  enthalpies  calculated  within 
the  local  density  approximation,  and  preliminary  predictions  of  the  dominant  ionization 
states  are  made.  Temperature-dependent  defect  formation  entropies  and  the  temperature 
dependence  of  the  preexponentials  are  incorporated  into  the  calculation  of  the  defect  den¬ 
sities.  Degenerate  Fermi-Dirac  statistics  are  used  for  the  electronic  equilibration,  and  the 
intrinsic  reaction  constant  as  a  function  of  composition  and  temperature  is  calculated. 
We  theoretically  substantiate  the  doubly  ionized  mercury  vacancy  as  the  dominant  defect 
in  HgCdTe,  and  expect  the  doubly  ionized  mercury  vacancy  densities  to  be  comparable 
in  HgZnTe.  We  predict  that  tellurium  antisites  are  donors  and  will  be  present  for  some 
annealing  conditions  in  sufficient  quantities  to  be  measured  and  possibly  to  affect  device 
performance. 


A-1 


1.  INTRODUCTION 


The  importance  of  native  defects  in  HgCdTe  is  undisputed,  with  the  dominant  defect 
believed  to  be  the  double  acceptor  Hg  vacancy.^  The  evidence  for  these  defects  is  largely 
indirect  and  depends  on  their  being  ionized  for  observation.  Neutral  defects  and  com¬ 
pensating  defects  are  more  difficult  to  measure,  and  neither  their  densities  nor  even  their 
presence  is  well  established. 

Our  goal  is  to  identify  the  important  native  defects  in  HgCdTe  as  a  function  of  temper¬ 
ature  and  mercury  pressure.  Predictions  of  the  absolute  defect  concentrations  are  difficult 
because  of  the  accuracy  required  for  reaction  enthalpies  and  entropies  that  enter  in  ex¬ 
ponentials.  Our  calculations  although  using  a  state-of-the-art  method,  are  subject  to  a 
number  of  limitations,  such  as  the  supercell  approximation  that  is  used  to  calculate  defect 
formation  enthalpies.  Despite  these  limitations,  for  which  error  bars  can  be  estimated, 
we  find  good  agreement  with  experiment  for  the  defect  densities  in  narrow-gap  HgCdTe. 
While  there  are  some  fitted  parameters  in  the  theory  —  for  example,  the  temperature 
variation  of  the  band  gap  —  none  have  been  chosen  to  fit  the  measured  defect  densities  we 
are  predicting.  Thus,  deviations  must  be  ascribed  to  deficiencies  in  our  approximations  or 
to  physical  mechanisms  that  have  not  yet  been  incorporated. 

In  our  previous  work,  we  used  the  linearized  muffin-tin  orbital  (LMTO)  method  within 
the  atomic  spheres  approximation  (ASA),’  and  later  the  full-potential  (FP)  Harris  Foulkes 
approximation.’  The  ASA  substitutes  a  spheridized  density  for  the  true  Honenberg-Kohn 
density  functional  in  the  local  density  approximation.  While  the  ASA  is  computationally 
fast,  it  cannot  reliably  predict  atomic  forces,  and  therefore  lattice  relaxations.  Because 
of  the  elimination  of  the  ASA  shape  approximation,  in  the  FP  calculations  we  have  pre¬ 
dicted  the  breathing  mode  relaxations  about  the  defect  sites.  Using  the  resulting  defect 
formation  enthalpies,  in  Ref.  3  we  predicted  ratios  of  neutral  native  defects.  In  the 
present  paper  we  have  made  several  important  improvements  to  this  previous  work;  (1) 
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a  gradient  correction  to  the  local  density  functional^  has  been  included,  which  we  find 
has  a  significant  impact  on  the  defect  formation  enthalpies;  (2)  absolute  defect  densities 
are  predicted;  (3)  the  reaction  constant  for  electron-hole  pair  production  as  a  function  of 
temperature  and  cadmium  composition  for  degenerate  Fermi-Dirac  statistics  is  predicted; 
and  (4)  the  primary  ionization  states  of  the  native  defects  have  been  tentatively  identified 
and  incorporated  into  the  defect  equilibrium. 

The  problem  of  predicting  the  defect  concentrations  is  complicated  by  the  fact  that  we 
are  dealing  with  an  alloy.  The  defect  formation  energies  were  calculated  for  HgTe,  the  pri¬ 
mary  component  of  the  pseudobinary  alloys  HgCdTe  and  HgZnTe.  Using  a  tight-binding 
model  we  have  shown‘  that  the  vacancy  formation  energy  has  a  nonlineu  dependence  on 
the  constituents  in  the  near  alloy  environment,  with  the  nonlinearity  being  most  dramatic 
for  the  removal  of  the  cormnon-lattice  atom  (i.e.,  tellurium  for  HgCdTe  and  HgZnTe). 
The  variation  of  the  mercury  vacancy  formation  energy,  while  less  dramatic  because  the 
first-nearest  neighbors  are  always  tellurium,  is  as  much  as  <^0.1  eV  (depending  on  the  con¬ 
stituents  in  the  second-neighbor  shell).  We  have  not  yet  incorporated  this  level  of  detail 
into  our  first-principles  calculations.  Additionally,  we  have  not  included  cadmium-  or  zinc- 
based  defects  (such  as  the  cadmium  antisite);  because  HgTe  is  the  dominant  constituent 
for  compositions  of  technologial  importance,  we  believe  this  exclusion  is  justified.  The  al¬ 
loy  effects  have  been  included  in  the  present  calculations  via  the  composition  dependence 
of  the  bamdgap  which,  because  of  its  temperature  dependence,  has  been  extracted  from 
experiment  and  the  shape  of  the  conduction  and  valence  bands. 

II.  REACTION  ENTHALPIES 

We  consider  the  following  native  defects:  mercury  vacancy  Vng,  tellurium  vacancy 
Vt*,  mercury  antisite  Hgj„  tellurium  antisite  Tewg,  mercury  interstitial  Hg;,  and  tellurium 
interstitial  Te/.  A  couple  of  corrections  have  been  added  to  the  reaction  enthalpies  calcu- 
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lated  within  the  H^l^^is-Foulkes  approximation  to  the  FP-LMTO,  First,  the  ASA  is  used 
to  determine  the  ionization  state  of  the  defect  —  that  is,  whether  it  is  a  donor  or  accep¬ 
tor.  We  have  not  yet  determined  the  ionization  energies  of  the  defect,  but  we  do  indicate 
whether  the  state  appears  to  be  “deep”  or  “shadlow.”  We  assume  that  only  the  shallow 
states  are  electrically  active.  Because  we  have  not  yet  determined  the  ionization  energies, 
we  will  for  the  present  assume  that  they  are  zero  for  the  shallow  states  —  that  is,  that  the 
donor  and  acceptor  levels  lie  close  to  the  conduction  and  valence  band  edges,  respectively. 
We  assume  that  the  deep  states  are  not  electrically  active.  This  approximation  is  likely  to 
be  good  because,  at  high  temperatures  where  the  defect  concentrations  are  equilibrated, 
carrier  concentrations  are  high  enough  so  that  free  carrier  screening  may  effectively  re¬ 
duce  the  shallow-state  activation  energies.  Details  of  the  calculational  method  used  to 
determine  the  ionization  states  of  the  various  defects  will  be  reported  elsewhere. 

The  second  correction  to  results  involves  the  incorporation  of  a  gradient  correction 
to  the  local  density  functional.  This  correction  has  not  yet  been  implemented  in  the  FP- 
LMTO,  and  thus  was  done  within  the  ASA.  We  expect  that  the  FP  gradient  correction  will 
be  nearly  equal  to  that  from  the  ASA,  because  the  density  gradient  is  predominantly  radial, 
the  nonspherical  components  eliminated  in  the  ASA  being  small.  In  Table  I  we  summarize 
the  formation  energies  including  the  gradient  correction  for  the  neutral  defect  reactions  in 
Eqs.  (l)-(6)  below.  Also  shown  are  the  cohesive  energies  of  the  constituent  compounds 
HgTe,  CdTe,  and  ZnTe  from  a  gradient- corrected  self-consistent  FP  calculation;  one  can 
see  that  the  agreement  with  experiments  is  quite  good.  The  tentative  identification  of 
defect  ionization  states  is  also  given  in  Table  I. 

III.  DEFECT  CONCENTRATIONS 

Defect  concentrations  are  determined  using  a  quasi-chemical  analysis  of  the  defect 
formation  reactions.*  We  consider  the  following  defect  reactions 


A-4 


HgTe-VH,Te  +  Hg, 

(1) 

Hg,  -  HgW. 

(2) 

2Hg,  -  HgHgx. 

(3) 

2HgTe  — ►  TeHgTe  +  2Hg^ 

(4) 

Hg,  Hg, 

(5) 

HgTe-Te,  +  Hg, 

(6) 

where  we  have  chosen  the  HgXe  unit  cell  and  &ee  atomic  mercury,  Hg,,  as  the  reference 
states.  Other  reactions  of  interest  can  be  obtained  by  taking  linear  combinations  of  these 
equations.  For  example  the  neutral  Schottky  defect  reaction  is  obtained  by  adding  Eqs. 
(1)  and  (2)  to  obtain 


HgTe-.VH,Te  +  HgVT.  (7) 

or,  as  it  is  more  commonly  written, 

0-Vh,  +  Vt..  (8) 

The  concentrations  of  the  defects  in  Eqs.  (l)-(6)  can  be  obtained  from  the  evaluation 
of  the  corresponding  reaction  constants: 


(9) 

=  [VV«]/PHg 

(10) 

=  [Hgx.]/p?., 

(11) 

=  [TeHgjpHg 

(12) 
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Kng,  =  [Hg/]/pHg 


(13) 


Krt,  =  [Te/]pHg-  (14) 

In  these  expressions,  pHg  is  the  mercury  pressure  in  atmospheres  and  square  brackets 
refer  to  concentrations  per  cubic  centimeter.  The  evaluation  of  these  reaction  constants 
was  discussed  in  our  previous  paper^  as  well  as  in  many  standard  texts. ^  In  the  present 
calculations,  we  have  used  the  formation  energies  from  Table  I  and  temperature-dependent 
entropies  given  in  Ref.  3. 

Ionized  defect  concentrations  can  be  determined  from  the  concentration  of  neutral 
defects  from 


[D  ]  =  {9D  l9D-)txp{iED-  -  Er)lkBT)  m  (15) 

for  donors  and 


[A']  =  {9a'/9a-)  exp((£r  -  E^.)/kBT)  [A‘]  (16) 

for  acceptors.  The  dot  and  prime  superscripts  correspond  to  a  positively  and  negatively 
charged  species,  respectively;  9x  is  the  degeneracy  of  the  state  X’,  Ed  and  Ea'  are  the 
positions  of  the  first  ionization  levels  for  the  donor  and  acceptor,  respectively;  in  the 
one-electron  picture,  Er  is  the  Fermi  energy,  iba  is  Boltzmann’s  constant,  and  T  is  the 
temperature  in  Kelvin.  Similar  expressions  are  obtained  for  the  second  ionization  state. 

In  addition  to  the  above  equations  for  the  determination  of  the  native  defect  popula¬ 
tions,  we  have  the  reaction  for  the  generation  of  electron-hole  pairs 

O-e'-fA.  (17) 

The  corresponding  reaction  constant  is 
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K,n  =  [h][e']=pn 


(18) 


where  p  =  [h]  and  n  =  [e'j,  as  in  the  usual  notation. 

In  general,  Kpn  depends  on  the  structure  of  the  conduction  and  valence  bands,  the  band 
gap  energy,  the  Fermi  energy,  and  the  temperature  variation  of  these  quantities.  Several 
limits  are  often  encountered  in  the  evaluation  of  Kpn-  First,  when  the  conduction  md 
valence  bands  are  parabolic,  although  not  necessarily  isotropic,  E  and  the  reaction 
constant  can  be  written  in  terms  of  the  Fermi-Diraw:  integrals  as 

^  (^)  (^)  (1®) 
where  Fj/j  is  the  Fermi-Dirac  function;  Ee,  and  Ep  are  the  conduction-band,  valence- 
band,  and  Fermi  energy,  respectively;  n.,.  and  m,  are  the  hole  and  electron  density-of-states 
effective  masses,  respectively;  and  h  is  Planck’s  constant.  In  the  nondegenerate  limit,  this 
reduces  to  the  familiar  expression 

if,,  =  *  (20) 

which  is  independent  of  the  Fermi  energy.  For  the  general  degenerate  case,  which  will 
apply  even  at  moderate  temperatures  for  narrow-gap  HgCdTe,  Kpn  will  depend  on  Ef 
and  therefore  in  general  will  not  be  independent  of  the  presence  of  extrinsic  carriers. 
Additionally,  the  assumption  of  parabolic  bands  may  be  poor  for  the  narrow-gap  materials, 
where  the  dispersion  near  the  conduction  band  edge  is  rather  more  linear  than  parabolic,* 
i.e.,  E  =  ak.  This  case  obtains  with  the  Fermi-Dirac  integral  function  of  order  2  and 


_^(2i:kBTmH\^'^  2  (kBT\^  rr  (Ef-E,\  p  ( 

-  2  p— j  F  [—)  [ 


Ev-  Ef' 


In  the  present  program,  we  have  used  Eq.  (21)  to  evaluate  Kpn-,  with  m*  =  0.43, 
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E,{x,  T)  =  Ec-E,=-  0.313  +  1.787x  +  0.444x*  -  1.237x»  +  0.932x" 


+  (0.667  -  1.714X  +  0.760x’)T/1000  (22) 

taken  from  Ref.  9,  and  a  chosen  to  yield  good  agreement  with  experimental  values  of  the 
intrinsic  carrier  concentrations.^® 

The  calculation  of  the  intrinsic  reaction  constant  for  narrow-gap  HgZnTe  is  more  dif¬ 
ficult  to  evaluate,  given  the  more  limited  data  base  for  evaluation  of  the  high-temperature 
band  gap,  effective  masses,  and  intrinsic  carrier  concentration.  From  the  zero-temperature 
band  structures,  the  valence-band  effective  masses  are  found  to  be  comparable  for  HgCdTe 
and  HgZnTe  with  equal  band  gaps,  while  the  conduction-band  effective  mass  for  HgZnTe 
is  slightly  larger  than  for  HgCdTe.^^  Overall,  though,  we  expect  Kpn  for  the  two  materials 
to  be  comparable. 

To  calculate  the  native  defect  concentrations  at  a  given  temperature  and  mercury 
pressure,  we  have  to  determine  the  Fermi  energy  which  satisfies  the  neutrality  condition 

2  riA'l  +  EWl  +  i*  1  =  2E  W'l  +  EWl  +  l''l-  (23) 

t  i  %  I 

The  concentrations  of  neutral  defects  which  are  independent  of  the  Fermi  energy  can  be 
solved  for  directly  for  a  given  png,  and  T. 

IV.  RESULTS 

The  low-temperature  (77  K)  hole  concentrations  as  a  function  of  mercury  pressure  for 
various  equilibration  temperatures  are  shown  in  Figure  1  for  x  =  0.2  Hg,Cdi_,Te.  Shown 
for  comparison  are  the  experimental  data  from  Vydyanath.^  In  the  calculations,  we  allowed 
both  atomic  and  electronic  equilibrations  at  the  high  temperature  at  which  the  annealing 
takes  place;  we  then  assume  the  total  defect  concentrations  are  frozen-in  upon  quenching. 


and  that  at  77  K  only  the  electronic  equilibrium  is  reestablished.  In  agreement  with 
Vydyanath,^  we  conclude  that  the  mercury  vacancy  is  the  dominant  defect  in  HgCdTe, 
and  it  behaves  as  a  double  accceptor.  We  have  not  yet  determined  the  activation  energy 
for  the  vacancy.  We  hnd  the  tellurium  antisite,  which  is  a  donor,  is  also  an  important 
defect.  In  Figure  2  we  show  a  breakdown  of  the  concentrations  of  various  defects  at 
500*  C  as  a  function  of  mercury  pressure.  At  low  mercury  pressures,  the  tellurium  antisite 
density  becomes  comparable  to  that  of  the  mercury  vacancy  and  results  in  a  high  degree 
of  compensation;  it  is  the  compensating  tellurium  antisite  donors  that  are  responsible  for 
the  roll-off  on  the  low-pressure  side  of  the  concentration  curves  in  Fig.  1. 

The  remaining  discrepancies  between  our  results  and  experiment  can  be  attributed  to  a 
number  of  factors.  First,  we  have  calculated  the  formation  energies  for  HgTe  and  applied 
them  directly  to  HgCdTe,  with  the  alloy  taken  into  account  in  the  calculation  of  the 
reaction  constant  for  electron-hole  pairs,  and  in  the  assumption  that  the  defect  ionization 
levels  are  at  the  band  edges.  Although  this  should  be  a  reasonable  approximation,  we  have 
shown  that  the  vacancy  formation  energies  are  sensitive  to  the  near-neighbor  environment.® 
The  removal  of  this  approximation  will  lead  to  a  raising  of  the  vacancy  formation  energy 
by  as  much  as  0.1  eV  and,  from  examination  of  Fig.  1,  will  yield  better  agreement 
with  experiment.  Next,  we  have  not  taken  into  account  the  nonradial  relaxation  about 
the  vacancy  site  and  its  effect  on  both  the  formation  energy  and  entropy,  nor  have  we  yet 
included  the  effect  of  the  entropy  of  local- mode  softening  of  the  ionized  vacancy  relative  to 
the  neutral  vacancy.  Additionally,  the  correct  ionization  energy  must  be  incorporated  into 
the  calculation.  Finally,  we  have  used  a  band  structure  in  calculating  the  reaction  constant 
for  electron-hole  pairs  which  yields  agreement  with  the  intrinsic  carrier  concentrations  and 
band  gap  at  T  <  400  K;  the  band  structure  at  higher  temperatures  is  speculative,  and  we 
need  further  experiments  or  theory  to  confirm  its  validity. 

At  the  annealing  temperature,  the  material  is  intrinsic  and  all  of  the  native  defects  are 


nearly  completely  ionized.  While  the  total  numbers  of  defects  sue  frozen-in  upon  quenching 
from  high  temperature,  the  ratio  of  ionized  to  nonionized  defects  does  change.  At  77  K 
the  material  is  extrinsic  and,  for  material  annealed  at  high  temperatures  (T  >  400* C),  is 
dominated  by  the  acceptor  level  of  the  mercury  vacancy.  Consequently,  upon  quenching, 
the  ratio  of  ionized  to  nonionized  defects  decreases  for  acceptor  defects  and  increases  for 
donor  defects. 

Defect  densities  present  after  a  low  temperature  (-  250* C)  annealing  under  mercury- 
saturated  conditions  are  of  technological  interest.  Otir  predictions  for  defect  concentration 
as  a  function  of  1  /T  for  mercury  pressures  corresponding  to  mercury  saturated  conditions 
are  shown  in  Figure  3.  At  T<  300®C,  a  mercury-staturated  annealing  step  is  effective  in 
reducing  the  vacancy  concentration  to  less  than  10^*cm"*.  Because  the  density  of  tellurium 
antisites  varies  as  compared  to  Phj  for  the  mercury  vacancy,  the  antisites  are  most 
important  on  the  tellurium- rich  side  of  the  stability  region,  and  thus  are  less  important 
for  these  mercury-saturated  annealing. 

While  the  mercury  vacancy  and  tellurium  antisite  are  the  primary  defects  that  will 
affect  the  electrical  activity  in  HgojCdo.jTe,  several  other  defects  are  of  interest  because  of 
their  role  in  diffusion,  for  example  the  mercury  and  tellurium  interstitials.  In  Figure  2,  we 
see  that  the  densities  of  mercury  antisites  and  mercury  interstitials  are  less  than  10^“cm"S 
while  the  densities  of  tellurium  vacancies  and  interstitials  are  completely  negligible.  The 
defect  formation  energies  used  for  the  calculation  of  the  interstitial  densities  were  from  our 
ASA  calculations  without  lattice  relaxations,  and  therefore  carry  a  larger  uncertainty  than 
the  antisite  and  vacancy  formation  energies.  We  expect  that  improvement  of  the  defect 
total-energy  calculations  will  result  in  a  potentially  significant  lowering  of  the  formation 
energy,  and  consequently  an  increase  in  the  interstitial  defect  concentrations.  Even  so,  it 
is  doubtful  that  we  will  find  tellurium  interstitials  at  densities  large  enough  to  account  for 
the  tellurium  diffusion,  which  was  believed  to  be  via  the  tellurium  interstitials. 


A-10 


While  the  tellurium  antisites  are  donors  and  are  predicted  to  be  present  in  significant 
densities  in  HgCdTe,  we  do  not  predict  that  equilibrated  tellurium  antisites  are  responsible 
for  the  experimentally  observed  p-to-n  conversion.  Even  if  such  a  p-to-n  conversion  were 
to  occur  as  a  consequence  of  refinements  in  our  calculations,  the  conversion  would  occur 
on  the  tellurium-rich  side  of  the  stability  region.  Because  antisites  are  likely  to  be  immo¬ 
bile,  relatively  high  densities  of  tellurium  antisites  introduced  at  high  temperatures  during 
growth  may  persist  through  low-temperature  annealings  undertaken  to  remove  mercury 
vacancies.  In  this  case  the  tellurium  antisites  could  be  responsible  for  the  p-to-n  conver¬ 
sion.  A  careful  correlation  of  the  n-type  character  with  the  high-temperature  preparation 
conditions  would  help  to  confirm  or  negate  this  possibility. 

Because  of  the  similar  magnitude  for  the  electron  and  hole  effective  masses  in  narrow- 
gap  HgCdTe  and  HgZnTe,  we  expect  that  the  intrinsic  reaction  constants  will  be  similar 
in  the  two  materials.  In  addition,  the  mercury  vacancy  formation  energy  will  be  similar 
in  the  two  materials  for  compositions  low  in  cadmium  or  zinc,  and  therefore  we  expect 
the  defect  concentrations  to  be  similar  in  the  low-i  materials  w’ith  comparable  bandgaps. 
The  tellurium  antisite  was  found  to  be  strained  in  the  HgTe  lattice,’-*  and  produced  a 
compressive  strain  in  the  surrounding  lattice.  Because  of  the  shorter  bond  lengths  in  the 
HgZnTe  alloy,  the  strain  associated  with  the  tellurium  antisite  is  larger  than  in  HgCdTe 
with  a  comparable  bandgap;  consequently  the  formation  energy  is  expected  to  be  larger 
in  HgZnTe.  Thus,  we  expect  the  tellurium  antisites  to  be  less  important  in  HgZnTe  than 
in  HgCdTe.  If  the  tellurium  antisite  is  responsible  for  the  residual  n-type  character  in 
low-temperature  annealed  material,  lower  n-type  carrier  concentrations  will  be  realizable 
in  HgZnTe. 

V.  CONCLUSIONS 

We  have  calculated  the  low-temperature  hole  concentration  as  a  function  of  anneal- 


A-11 


ing  temperature  and  pressure  for  HgCdTe.  We  find  good  agreement  with  the  fit  to  the 
experimental  data.  Remaining  discrepancies  can  be  attributed  to  relaxation  and  finite 
ionization  energies  in  the  vacancy  formation  energies,  and  uncertainties  in  Kp„  at  high 
temperatures.  We  conclude,  in  agreement  with  experiment,  that  the  dominant  defect  is 
the  doubly  ionized  mercury  vacancy.  We  do  not  make  any  conclusions  about  the  pres¬ 
ence  of  the  singly  ionized  Hg  vacancy;  we  will  examine  the  activation  energies  for  both  the 
singly  and  doubly  ionized  states  in  future  work.  We  find  that  tellurium  antisites  are  donors 
and  Me  present  in  significant  densities  for  low  mercury  pressures  at  temperatu'es  above 
350  K.  At  lower  temperatures  and  high  mercury  pressures,  at  which  annealing  typically 
takes  place,  tellurium  antisite  equilibrium  densities  are  negligible.  However  if  their  den¬ 
sity  remains  at  levels  corresponding  to  high  temperature  because  their  low  temperature 
diffusivities  are  low,  then  such  antisites  could  affect  mobilities  and  even  be  the  “residual 
donor."  None  of  the  other  native  defects  are  present  in  sufficient  quantity  to  influence 
device  characteristics. 
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FIGURE  CAPTIONS 


FIGURE  1.  Hole  concentrations  at  77  K  as  a  function  of  mercury  pressure  for  material 
annealed  at  high  temperatures.  Experimental  data  were  taken  from  Vydyanath^ 
Theory  is  shown  as  soHd  lines  for  pressures  within  the  stability  region  at  a  given 
temperature. 

FIGURE  2.  Predicted  defect  densities  at  500  K  as  a  function  of  mercury  pressure 

FIGURE  3.  Defect  densities  as  a  function  of  temperature  for  mercury-saturated  an¬ 
nealing  conditions 
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TABLE  I.  Formation  energies  and  ionization  states  of  the  native  point  defects.  Formation 
energies  refer  to  the  neutral  defect  reactions  in  Eqs.  (l)-(6)  in  text.  See  text  for  discussion 
of  shallow  and  deep.  Experimental  cohesive  energies  are  shown  in  parentheses. 


Defect  Formation  energy  (eV) 

Ionization  state 

Viig 

1.93 

shallow  acceptor 

Vt. 

2.39 

donor 

TeHg 

2.68 

shallow  donor 

HgT. 

0.75 

deep  acceptor 

Te; 

4.47 

shallow  donor 

Hg; 

1.75 

shallow  donor 

ZnTe 

4.66  (4.8) 

CdTe 

4.17  (4.4) 

HgTe 

3.37  (3.3) 
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DISLOCATION  SPACE  CHARGE 


DISLOCATION  SPACE  CHARGE 


Represent  the  dislocation  core  by  a  cylindrical  charge  density  -p  in  a  cylinder  of  radius  a 
in  n-type  material  with  donor  density  n.  The  potential  distribution  in  the  depletion  approxima¬ 
tion  is  given  by  the  Poisson  equation  with  charge  density  p 

p  =  (.p  +  n)e  r  <  s 

=  ne  a  <  r  <  R 

=  0  R<r 

where  R  is  the  depletion  radius,  related  to  a  by  the  charge  neutrality  requirement  pa^  =  nR^.  By 
integrating  the  Poisson  equation  we  find  for  the  potential 

<|)(r)  =  <t)(a)  +  g  [^R^ln  J  ^  {r^  -  a^)] 

and  we  take  the  potential  in  the  bulk  (r  >  R)  as  reference;  we  can  evaluate  the  potential  at  the 
core  radius 

<t>(a)  =  ^  [2  • 


In  MCT  a  reasonable  value  of  the  dislocation  core  radius  is  5  x  10'*  cm  and  a  potential  relative  to 
the  bulk  at  the  core  radius  of  0.1  V.  A  table  of  depletion  radius  versus  doping  level  for  these 
values  is 


n,  cm'® 

10’® 

3x  10’® 

ioi4 

3x10’^ 

10’5 

3x10’® 

10’® 

3x  10’® 

R.  cm 

I.TOx  10'^ 

1.01  xKT^ 

RMI9 

mm 

4.32  X  10-® 

These  radii  are  close  enough  to  the  Debye  length  to  make  the  depletion  approximation  question¬ 
able.  The  exact  space  charge  equation  is 


iiLr^  = 
r  dr  dr 


[sinh((p  -  (Pf)  =  sinhtpp] 
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where 


f 


9 


9 


<P  =  S  L"  =  -PF  =  ^  =  sinh-' 


kT 


e^n, 


kT 


For  <p  small  this  is 


I  Ll 

with  the  extrinsic  Debye  length  Ln  =  — ^ ^ 

cosh9F  e^ND  Np 


,  Nd  in  cm-3 


Then  with  x  =  we  have  x  ^  =  9,  the  zero-order  modified  Bessel  equation.  The 
Ld  X  dx  dx 

solution  of  interest  is  the  decreasing  Hankel  function  Ko(x)- 


The  depletion  approximation  should  be  valid  for  r  <  R  =  Lp.  This  does  not  leave  much 
room  since  R  and  Ld  are  comparable.  One  way  to  proceed  is  to  let  9  =  CKo(x)  and  to  evaluate  C 
from  9(a)  at  x  =  a/Lo-  Then  we  can  determine  the  effective  scattering  radius  of  the  dislocation  as 
the  radius  at  which  the  field  is  reduced  to  a  small  enough  value,  say  10*2  V/cm.  The  space 
charge  field  is  given  by 


d9  1^_^_  kT  d9 
dr  e  dr  eLp  dx 


\J  ^  CKi(x)  =  4.75  X  10-5  ^  CK,(x). 


For  the  core  radius  of  5  x  10*^  cm  and  core  edge  potential  of  0. 1  V,  and  a  donor  concen¬ 
tration  of  lO^'*  cm*3  we  find  Ld  =  5.26  x  10*5  cm,  a/Lo  =  9.5  x  10*^. 

For  X  small,  Ko(x)  =  In  —  =  0. 1 16  -  Inx,  Kq  [7^]  =  0. 1 16  +  6.96  =  . 

yx  \^d)  C 

Since  <j>(a)  =  ^  9(a)  =  0.0259(a)  =  0. 1  this  gives  C  =  0.56,  and 

C 


E  =  4.75  X  10'^  X  10^  X  0.56  K,(x)  =  266  Ki(x) . 


For  X  large,  E  =  266>^/  e  *  =  333  -7= 

2x 


To  reduce  E  to  10*2  V/cm  we  need  e’VV^  =  3  x  10’^,  x  =  9.3,  r  =  4.9  x  10“^  cm. 
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PIEZOELECTRIC  FIELD  OF  A  STRAIGHT  DISLOCATION 
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PIEZOELECTRIC  FIELD  OF  A  STRAIGHT  DISLOCATION 


The  polarization  produced  by  a  strain  ujk  is  given  by  Pi  =  -eijkUjk,  where  eyk  is  the  ijk 

component  of  the  piezoelectric  tensor  (piezoelectric  coupling  constant)  and  where 
1  fduj  dujA 

Ujk  =  —  I +  -^1 ,  with  uk  the  k-component  of  the  stress-induced  distortion. 

In  the  sphalerite  structure  the  only  nonzero  component  of  the  piezoelectric  tensor  is  ei4  if  i,  j,  and 
k  are  all  different.  Since  there  is  no  free  charge  in  an  insulating  crystal,  the  displacement  is 
solenoidal,  divD  =  0,  and  the  potential  is  given  by  V2<)>  =  -  4jtp(r),  with  an  effective  bound 

charge  density  p(r)  =  ^ — [ejjkUjkCr)] . 

£  »Vv. 


If  we  take  the  z-axis  along  the  direction  of  a  straight  dislocation,  all  the  strains  are  functions  of  x 
and  y  only.  For  an  edge  dislocation*  the  only  shear  strain  is 

u,  bx(y2  -  x2) 

^  471  (1 -V)(x2  +  y2)2 

where  b  is  Burger’s  vector  and  v  is  Poisson’s  ratio.  For  a  screw  dislocation  along  the  z-axis  with 
b  =  bz  the  shear  strains  are 

u  u  bx 

^  47C  (x2  +  y2)  ’  4%  (x2  +  y2) 

It  is  apparent  that  no  bound  charge  density,  and  hence  no  piezoelectric  field,  is  associated  with  a 
pure  edge  dislocation.  The  bound  charge  density  associated  with  a  pure  screw  dislocation  is 

p-beu/^  x  ^  y  ]_bei4  y^-x^ 

\dx  +  y^  dy  x2  +  y2  j  27t  1^2  +  y2)2 


*  See,  ios  example,  D.  Hull  and  D  J.  Bacon,  “Inuoduc^on  to  Dislocations,”  3rd  ed.  Pergamon  1984. 


